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Using flow visualizations in Couette geometry, we demonstrate the existence of Taylor-like vortices 
in the shear-banding flow of a giant micelles system. We show that vortices stacked along the 
vorticity direction develop concomitantly with interfacial undulations. These cellular structures are 
mainly localized in the induced band and their dynamics is fully correlated to that of the interface. 
As the control parameter increases, we observe a transition from a steady vortex flow to a state 
where pairs of vortices are continuoulsy created and destroyed. Normal stress effects are discussed 
as potential mechanisms driving the three-dimensional flow. 

PACS numbers: 47.50.-d, 47.20.-k,83.60.-a, 05.45.-a, 83.85.-Ei 



In contrast to Newtonian fluids, complex fluids are 
structured at the mesoscopic scale, typically from 10 nm 
to 100 /im, corresponding, for example, to the mesh size 
in entangled polymers solutions or to the bubble size in 
foams. Such a length scale provides internal degrees of 
freedom with characteristic times in the order of 10~^ to 
100 s (position, orientation, composition, ...) that can 
easily be excited by the flow, often resulting in a new or- 
ganization of the structure of the fluid and of the flow it- 
self. The relation between stress and shear rate a = f(^) 
becomes strongly nonlinear and the flow evolves towards 
a spatially heterogeneous shear-banded state where re- 
gions of differing viscosities coexist. The shear-banding 
phenomenon is ubiquitous in complex fluids, from surfac- 
tants solutions m [2] , to telechelic polymers [3l [4] , emul- 
sions [5l [6] , granular materials [7 and foams [8] . Among 
these systems, surfactant giant micelles have generated 
an abundant literature [1] and are often considered as 
model systems for the study of shear-banded flows. The 
(J = /(7) curve is made of two stable shear-thinning 
branches, characteristic of two distinct states of the sys- 
tem, separated by a stress plateau (a = dp) that extends 
from 7^ to 7/j,. In the plateau regime, the flow is unstable 
and evolves towards a banded state where the viscous and 
fluid phases coexist at constant stress cFp [9] . This picture 
is supported, among others, by investigations focused on 
the determination of the flow field [IQl El H [E], re- 
vealing a discontinuity of the velocity profile in the flow 
plane. These ID measurements also allowed to estab- 
lish the existence of large spatio-temporal fluctuations of 
the local flow fleld [141 HSj^ suggesting the presence of 
three-dimensional flow [15]. Recently, using the scatter- 
ing properties of the induced band, we showed that the 
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interface between bands is not flat as usually assumed 
but undulates along the vorticity direction [16, 17 . Such 
a scenario has been qualitatively reproduced by a nu- 
merical perturbation analysis of the modifled Johnson- 
Segalman equation [18]. In addition, the model predicts 
the formation of velocity rolls stacked along the vorticity 
direction. 

In this letter, we report the flrst direct observation of 
Taylor-like vortices in the shear-banding flow of a semi- 
dilute micellar solution. We follow simultaneously the in- 
terface and flow dynamics and show that, as the interface 
between bands becomes unstable, a secondary axisym- 
metric vortex flow develops. The Taylor-like vortices are 
stacked along the vorticity axis and are mainly located 
in the high shear band, their size depending on the ap- 
plied shear rate. The radial velocities are found to be 
two orders of magnitude lower than the azimuthal base 
flow. Using flow visualizations, we show that the inter- 
face and flow dynamics are strongly correlated : as the 
control parameter increases, we observe a transition be- 
tween a steady vortex flow towards a state dominated by 
continuous creations and annihilations of cellular struc- 
tures. 

In our experiments, the working fluid is a 11% wt. semi- 
dilute aqueous mixture of surfactant (cetyltrimethylam- 
monium bromide) and salt (sodium nitrate) at 28 ° C. It 
is contained in a transparent Perspex cylindrical Cou- 
ette device (gap e=1.13 mm, inner radius i?i=13.33 
mm, height h=4Q mm) designed for imaging of the 
flow gradient-vorticity plane and side views of the cell 
(Fig. [l a). The device is attached to the shaft of a com- 
mercial rheometer (Physica MCR500), with the inner 
cylinder rotating. The experiments are performed with 
the shear rate as the control parameter and the response 
of the system is recorded from two distinct channels, in 
a simultaneous way. The flrst one is dedicated to flow 
vizualizations in the flow- vorticity (yiQ.,\iz) plane using 
seeding anisotropic reflective particles (anisotropic mica 
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FIG. 1: (a) Experimental setup for observations in the 
(uojUz) and (ur,Uz) planes, (b) Steady-state flow curve and 
corresponding banding structure in the (ur, u^) plane for dif- 
ferent shear rates along the stress plateau. The left and right 
sides of each picture corresponds respectively to the inner and 
outer cylinders. 



platelets from Merck at a volume fraction of 6.10"^), as 
it was done for the study of the dynamic regimes of iner- 
tial Taylor- Couette flows [19]. The fluid is illuminated by 
ambient light and the intensity Io{z) reflected in the ve- 
locity gradient direction is collected on a digital camera 
(CCDl). The second channel is dedicated to the obser- 
vation of a cross-section {{ur^Uz) plane), using either a 
He-Ne laser sheet (A=633 nm, thickness 1.5 mm) prop- 
agating along the velocity gradient axis or white light. 
A digital camera (CCD2) records the scattered inten- 
sity at 90°. This channel provides information on the 
interface dynamics and allows the computation of the 
flow field using particle image velocimetry (PIV). Figure 
[T]b shows the steady-state shear stress as a function of 
the shear rate for the sample sedded with mica flakes. 
This flow curve, measured in strain-controlled mode, ex- 
hibits a stress plateau with a significant positive slope 
at dp = 150 ± 1 Pa, extending from 7/ = 7 ± 0.5^"^ 
to = 95 ± 5s~^. Along the stress plateau, the sys- 
tem is organized into two bands separated by an inter- 
face that undulates along the vorticity direction with a 
well-defined wavelength [161 [III • The proportion of the 
induced band and the wavelength of the undulations are 
increasing functions of the applied shear rate (inset in 
Fig. [l]b). The flow structure associated with the interfa- 
cial undulations is illustrated in Fig. [2] for 7=40 s~^. On 
the snapshot of a cross-section taken when steady-state 
is achieved, we overlay the corresponding instantaneous 
flow field computed using PIV (see Fig.[2la). If, at first 
sight, no particular pattern can be extracted, axial and 
radial velocity components are detected, suggesting that 
the flow is not purely orthoradial. The existence of a sec- 
ondary flow is confirmed by side views of the Couette cell 
under ambient illumination that reveals a succession of 
axisymmetric dark and bright stripes stacked along the 
vorticity axis (Fig. [2jb). Such a picture is reminiscent of 
the primary pattern of the classical Taylor instability [19] 



and points out that the base flow is replaced by a vor- 
tex flow. In fact, the variations of the reflected intensity 
results from differences in the orientation of the seeding 
flow-aligned platelets : dark stripes are associated with 
regions where the flakes do not reflect the light along the 
observation direction, indicating radial flow, while bright 
stripes correspond to flow perpendicular to the observa- 
tion direction [20 . In other words, the vortices axis are 
centred on bright stripes, the radial flow coming from 
the adjacent dark stripes. The comparison between the 
reflected light distribution Io{z) and the interface profile 
noted I{z) (Fig. [2]d and c) shows that the flow is ori- 
ented radially at the crests level of the interface profile. 
This means that, the interfacial wavelength A is twice the 
wavelength A/^ characterizing the axial structure of the 
flow. To differentiate the inflow and outflow and quan- 
tify the corresponding velocities, we use the flow field 
infered from PIV. As mentioned above, due to experi- 
mental limitations, it is not possible to distinguish cellu- 
lar structures from our measured flow field. Nonetheless, 
to test its validity, we compute the following quantity 
that provides a measurement of the mean local orienta- 
tion of velocity vectors with respect to the vorticity di- 
rection : Ip{z) = (^\vz{z,ri)\ /y/v'^{z,ri) -^v'^{z,ri)y 

Taking into account the relation between the intensity 
profile Io{z) and the radial velocity component [21 , we 
expect the mean orientation profile Ip{z) to be similar to 
Io{z). Fig. |2j(d-e) show that the axial structure Ip{z) is 
consistent with the intensity profile Io{z) gathered from 
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FIG. 2: (a) Stationary banding structure at 7=40 s~^ and cor- 
responding velocity field computed from the PIV algorithm in 
the (r, z) plane, (b) Side view of the Couette cell under ambi- 
ent illumination, (c) Interface profile I(z) extracted from (a), 
(d) Intensity distribution Io{z) along the vertical axis inferred 
from (b). (e) Mean orientation Ip{z) of velocity vectors with 
respect to the vertical axis, (f) Radial velocity profile Vr at a 
position r = + 0.5e in the gap highlighting the inflow and 
outflow, (g) Direct observation in the plane (r, z) of vortices 
stacked along the vorticity direction. In this last case, the 
sample is free of tracers and illuminated with white light. 
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side views of the Couette ceh, indicating that, even if our 
computed instantaneous flow field is noisy, some relevant 
information can be extracted from the latter. Fig. [2]f 
displays the radial component Vr{z) at a given position 
close to the interface between bands. First, we find that 
the outflow and inflow are aligned respectively with the 
outer and inner crests of the interface proflle. Second, 
radial velocities are in the order of 0.1 mm.s~^, namely 
at least two orders of magnitude lower than the mean 
base azimuthal velocity. The weakness of the secondary 
flow makes the determination of the tracers displacement 
diflicult and explains the absence of pattern in our map- 
ping of the local velocities in the cross-section. 
From those experiments we can conclude that vortices 
develop in the shear-banding flow of giant micelles. How- 
ever, they do not allow the determination of the radial 
structure of the vortex flow. To address this question, we 
take advantage of the scattering contrast between both 
bands, due to concentration effects. Using white light il- 
lumination of the sample, free of seeding particles, we are 
able to directly image the cellular structures (Fig. |2]g). 
The latter are detectable because of local index variations 
probably due to the concentration difference between the 
two phases. The counter-rotating vortices form concomi- 
tantly with the interface undulations and are stacked 
along the vorticity direction. They grow from the in- 
flow regions and flnally span the high shear rate band 
with a reduced extension in the low shear region when 
the steady-state is achieved (see Movie in Suplemental). 

To investigate the role of the control parameter on the 
flow dynamics, we record the interface proflle I{z) and 
the reflected light distribution Io{z) as a function of time 
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FIG. 3: Evolution in gray scale of the interface profile I{z) and 
the reflected intensity distribution Io{z) as a function of time 
during step shear rate from rest to (a-b) 7=35 s~^ and (c-d) 
7=75 s~^. The black and white dots observable in Figs. 3.c 
and d correspond to a bubble, trapped in a recirculation zone 
and that periodically crosses the observation fields. 



and construct in gray scale the corresponding spatiotem- 
poral diagrams I{z, t) and /o(^, t) • Note that the /o(^, ^) 
diagrams (Fig. |3]b and d) are fully reproduced by plot- 
ting as a function of time, an intensity profile along a 
vertical line in the fiow cross section going through cel- 
lular structures in Figure 2.e. Fig. [3]a and b show the 
patterns observed for 7=35 s~^. After a transient of du- 
ration r corresponding to the formation of a fiat interface 
between bands, we observe the destabilization of the in- 
terface together with the development of axisymmetric 
cellular structures. Our experimental resolution do not 
allow to determine whether the Taylor-like velocity rolls 
are triggered by the interfacial instability or vice-versa. 
The growth of the interface instability and the appear- 
ance of radial and axial components in the fiow field are 
followed by saturation in a steady-state characterized by 
well-defined wavelength and amplitude of the interface, 
and a steady vortex structure. The patterns are only 
slightly modulated in time with a period T. This mod- 
ulation corresponds to a small variation of amplitude in 
the inner crests (dark regions in Figs, [sja and c)) of the 
interface profile [T^ and produces a rope-like structure in 
the spatiotemporal sequence of the refiected light Io{z). 
It could result from periodic oscillations of the vortices 
boundaries. These extremely regular patterns are rep- 
resentative of the system dynamics on a large part of 
the stress plateau, for 7 typically between 15±1 s~^ and 
60±5 s~^. In this range of 7, we find that the asymptotic 
wavelength of the patterns increases with the shear rate 
(Fig.[4ja) while the period of the modulation keeps a con- 
stant value (Fig.|4]b). For comparison, the wavelength of 
the spatial structure Ip{z) obtained from PIV measure- 
ments has been added in Fig.|4]a. For all shear rates, the 
height of a vortex pair and the interfacial wavelength are 
in good agreement. The circulation cells essentially span 
the high shear rate band in the radial direction, but their 
extension along the vorticity direction can reach about 
two gap widths as the control parameter is increased. 
For higher values of the shear rate (7 > 65 s~^), the sys- 
tem exhibits a more complex dynamics as illustrated in 
Fig.|3]c and d. Once again, the dynamics of the tracers is 
remarkably correlated with the one of the interface. Af- 
ter the transient, the interface amplitude saturates but 
the wavelength continuously evolves with time due to nu- 
cleation and merging of inner crests. Since a pair of vor- 
tices develops between two consecutive inner crests, such 
a dynamics corresponds to a succession of creations and 
annihilations of pairs of vortices. These processes do not 
tend toward a stationary situation, rather leading to a 
chaotic dynamics. Note that at low shear rates (below 
15 s~^), the reduced size of the induced band prevents 
the description of the fiow dynamics. 
The present results demonstrate that the shear-banding 
fiow of our semi-dilute giant micelles solution is three- 
dimensional, far away from the classical picture usually 
used in this field [TOl |TT] . Very recently, the formation of 
Taylor-like velocity rolls has been predicted by a nonlin- 
ear numerical study of the modified Johnson- Segalman 
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FIG. 4: (a) Comparison between the asymptotic interfacial 
wavelength A (□) and the height of a vortex pair inferred 
respectively from Io{z) (■) and Ip{z) (•). (b) Period T of the 
temporal modulation of the interface profile (o) and of the 
rope-like structure of the reflected intensity (■). The error 
bars represent the standard deviation obtained over several 
experiments. The time and length scales are computed by 
direct Fourier transform. 



model [T8^ : the 3D flow is triggered by the instability 
of the ID gradient banded state with respect to interfa- 
cial undulations along the vorticity direction, with nor- 
mal stress and shear rate jumps accross the interface as 
the underlying mechanism [22]. Some of our observa- 
tions are qualitatively captured by this model. How- 
ever, it does not reproduce the complex flow dynamics 
we report. Alternatively, another possible origin for the 
Taylor-like vortices is the elastic instability originally ob- 
served in dilute polymer solutions [23^ ^ and more re- 
cently in rod-like systems [25l [26] . In working conditions 
where the streamlines are curved and the fluid is elas- 
tic, the base orthoradial flow can become unstable due 
to a centripetal body force called 'hoop stress' [23 . Such 
an hypothesis has been proposed to explain fast tempo- 
ral fluctuations of the flow field that lead to suspicions 
for the presence of 3D shear-banded flow in a concen- 
trated giant micelles system [15 . Following the analysis 
in Ref. [15], we compute a modified Weissenberg num- 
ber S = {e/Ri)^^'^Wi as a function of the applied shear 



rate 7 . Here Wi is an approximate Weissenberg number 
obtained from the stress ratio Ni/a and (e/Ri) the 
gap ratio. The critical value above which an elastic 
instability appears is around 6 [23^, [28]. This stability 
criterion has been established in the framework of the 
Oldroyd-B model. If we assume that each band can be 
modelled by an Oldroyd-B equation [29 , we can compare 
the modified Weissenberg number S in each band with 
the critical value Note that, if the Oldroyd-B model 
does not describe shear-thinning properties, the shear- 
thinning character is supposed to reduce the instability 
threshold [30 . In the linear regime, S = and it is found 
to increase with 7 in the nonlinear regime, namely above 
1 s~^ (data not shown). Along the first branch, S never 
exceeds 0.5 while for 7 > 7/j,, S > Sc, suggesting that the 
induced band is likely to undergo an elastic instability, 
leading to formation of cellular structures. 
In summary, our results provide the first direct evidence 
for the existence of Taylor-like velocity rolls in the shear- 
banding flow of wormlike micelles. Moreover, we showed 
that the dynamics of the Taylor-like vortices is strongly 
correlated with the interface dynamics, leading to funda- 
mental questions about the physical mechanism driving 
the three dimensional flow. Between the interface insta- 
bility and the formation of Taylor-like vortices, one chal- 
lenging task is now to determine which one of these phe- 
nomena drives the other, and under which physical mech- 
anism. More generally, our results could potentially ex- 
plain some fluctuating behaviors observed recently in gi- 
ant micelles systems and open new perspectives towards 
the understanding of shear-induced instabilities and het- 
erogeneous flows in complex fluids. 
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